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We analyze by a renormalization method, the dynamics of a particle in a infinite square-well 
potential driven by an external monochromatic field. This method set up for Hamiltonian systems 
with two degrees of freedom allows us to analyze precisely the stability of the trajectories of the 
particle as a function of the amplitude e of the external field. We compute numerical values of 
e for which the motion of the particle with frequency uj is broken and a transition to a chaotic 
behavior occurs. We obtain the critical function eduj) associated with this system as a function of 
' the parameters such as the frequency of the field and the width of the potential. 
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I. INTRODUCTION 



Due to the existence of as many conserved quantities as degrees of freedom, the trajectories of an integrable 
Hamiltonian system are confined to evolve on invariant rotational tori. On a given torus, the dynamics is regular, 
i.e. conjugate to a hnear flow with frequency w in action-angle coordinates. This regularity is broken by any small 
perturbation, and chaotic trajectories appear : The phase space of a Hamiltonian system close to integrable is in 
general a mixture of regular and chaotic motions. For Hamiltonian systems with two degrees of freedom, these 
Ph ' invariant rotational tori, also called KAM (Kolmogorov-Arnold-Moser) tori, act as barriers in phase space. 
— < For a given one parameter family of Hamiltonians {-ffe} with i?e=o integrable and e the amplitude of the perturbation, 
it appears from numerical evidences that for e smaller than a critical value denoted £c(w), there exists a KAM torus 
with frequency lo, and this torus is broken for larger values by resonance or overlapping of resonances. The function 
u) I— > ec(w) is called the critical function associated with the one-parameter family of Hamiltonians {-ffe}. This function 
^ , contains the information on the existence of invariant tori in phase space as the parameter e increases. In particular, 
for £ > supj^g^ ec(<^), there is no longer any KAM torus in the region of phase space corresponding to the set of 
J;0 ' frequencies B, and we have large scale stochasticity. 

' Renormalization methods have been defined and studied numerically for the analysis of stability of Hamiltonian 
I systems with two degrees of freedom ■ The aim is to describe the break-up of a given invariant torus. The idea 
^-H is to set up a transformation that focuses on a specific region of phase space around the given torus. It acts as a 
microscope in phase space, looking at the system on smaller scales in phase space and on longer time scales. The 
, complete renormalization method PJ|,^ is a tool to compute precisely the critical function £c{i^)- It has been verified 
• ^ that for specific models (like a forced pendulum) the critical couplings obtained by renormalization coincide with 
other methods like Greene's residue criterion or Laskar's frequency map analysis 

The system we analyze is a particle of mass m in an infinite square-well potential Vsq of width 2a driven by an 
external monochromatic field with amplitude £ and frequency Q. The Hamiltonian of this system with 1.5 degrees of 
freedom is the following : 



H{p,x,t) ^ f- + VsQix)+£xcos{nt), (1.1) 

where 

Vsq{x) — for |a;| < a and Vsq{x) — +co for > a, 
and £ > 0. This system has been studied in Refs. [01T2] by approximate renormalization methods. Without external 



field, i.e. for e = 0, the system is integrable and the motion is periodic with frequency oj — y 2mi? "^here E denotes 
the energy of the system. For e 7^ 0, some of these regular motions disappear. In particular, there are resonances 
when the frequency of the external field U, is commensurate with the frequency of the motion, i.e. when 17 is equal 
to (P and Q are relatively prime integers). The interaction of these resonances breaks up some invariant tori (in 
the spirit of Chirikov's criterion |13 1^). The critical function £c{i-o;rn,a,Vl) is the critical value of the amplitude of 



the field for which the motion with frequency uj is broken. 
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The aim of this paper is to apply renormaUzation to a specific model and to compute numerically Sc{oJ', m, a, We 
use this critical function to locate chaotic zones and to determine critical parameters for which large scale stochasticity 
occurs. We compare some of the results obtained by the renormalization transformation with the ones obtained by 
other existing methods such as Greene's residue criterion, in order to validate the resul ts o btained by renormalization. 

In Sec. ^ we explicit the model and give some information on its dynamics. In Sec. Ill, we give a short description 
of the renormalization method, and in Sec. IV we compute numerically and analyze the critical function ec('^; fn^ ol, 



II. MODEL 



The external field induces resonances in the system when 17 — ^lu where uj is the frequency of the unperturbed 
motion, Q the one of the external field and P, Q relatively prime integers. Due to the specific form of the interaction 
between the particle and the field, the largest resonances (i.e. of order of the amplitude e of the field) are obtained 



with Q = 1 and P odd. This can be seen by writing Hamiltonian (1.1) in action-angle variables |l/ 

H{A, ^, t) = -^A' _ 1^ ^ J_ cos(n^ - nt). (2.1) 

Ti odd 

There is a resonance when — fJ, which corresponds to nu!=Jl since ip — ^ is the frequency w of the motion. This 



resonance will be denoted l:n in what follows. Hamiltonian (2.1) can be mapped into a time-independent Hamiltonian 



with two degrees of freedom by considering —ilt as a, new angle variable 



H{Ai,A2,(pi,Lp2) = ^Al~nA2 2- V -^cos{n(pi+(p2). (2.2) 

n odd 

We rescale time by a factor i.e. we multiply Hamiltonia n (l2.2| ) by We notice that this rescaling of time changes 
the frequency of a quasiperiodic motion of Hamiltonian (|2.2D by a factor fi. The new rescaled frequency is now cu/fl. 
We rescale the action variables by replacing H{A,cp) by X~^H{A/X,(p) with a factor A = 7r^/(4ma^f2). We not ice 
that the rescaling in the actions does not change the equations of motion. After this rescaling, Hamiltonian ( |2.2| ) is 
equal to 

^2 2 

H{A, y,)^^-A2-e'}2^ cos(n(^i + (^2), (2.3) 

n odd 

where e' is the dimensionless amplitude of the external field given by : 

e 



(2.4) 



For e' = 0, Hamiltonian ( |2.3| ) depends only on A — {Ai,A2) and the equations of motion show that Ai(t) and 
A2{t) are constant, and fi{t) = ujt + ipi^o and ip2{t) — ~t + ip2.o- The trajectories of the system ( |2.3D evolve on 
three-dimensional energy surfaces in the four-dimensional phase space, and for e' = the trajectories with frequency 
UJ are confined to evolve on a two-dimensional torus (on the energy surface) with frequency vector u) = {lu,—1). 
For s' > 0, this system has an infinite number of main resonances (given by the condition nipi+ip2 ~ 0) located around 
Ai = l/n where n is an odd integer, accumulating at Ai — 0. The width of the n-th resonance zone is approximately 
equal to Ay/P/n [0. For large values of e' the torus with frequency vector lj is broken by overlapping of resonances. 
In order to have an estimate of the critical value of s' of the break-up of the invariant torus with frequency lo, we apply 



Chirikov's criterion |13|. For a torus with frequency lo located between the two primary resonances l:n and l:ri-|-2, 
the overlapping is obtained when the sum of the two half-width of these two resonances is equal to their distance, i.e. 
for 

= T?-^-^- (2-5) 

4(n+l)2 ^ ^ 

This value overestimates the critical values of the threshold of the break-up as it has been noticed in Ref. jl^ for 
this model. A convenient way to compute the value of e' for which the torus is broken is to set up a complete 
renormalization in the spirit of Refs. p],p^. 
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III. RENORMALIZATION METHOD 



In this section, we give a description of the renormahzation method we apply to the model described in the previous 
section. 

First, we shift the actions such that the invariant torus with frequency lo is located at Ai = for the unperturbed 
Hamiltonian (for e' = 0) : A'^ = Ai — to and A2 = A2. Hamiltonian ( p7^ ) is equal to : 

H{A',ip)=i^-A' + ^{n-Ar-e' Yl ;Jlcos(«..<^), (3.1) 

n odd 

where ft = (1, 0) and w = (w, —1). 

The renormahzation transformation is a map H' = Tl{H) acting within the family of Hamiltonians of the form : 

H{A,ip)=u;-A + V{n-A,ip), (3.2) 



and Hamiltonian (3.1) is an element of this family. Furthermore Hamiltonian (3.1) is quadratic in the actions; this 
feature is useful for a simplification of the implementation of the transformation (see Step 4 below) Q. However, a 
similar version of the transformation (obtained by slightly changing the way Step 4 is implemented) can be defined 
and studied numerically |^,|o|. 

The transformation TZ is based on the continued fraction expansion of lo : 



ao 



i r 

i — = ao, 0,1, 



ai+--- 

The best rational approximates of lo are given by the truncations of its continued fraction expansion : 

— = [ao,ai, ... ,afc = ooj. 
Qk 

The corresponding periodic orbits with frequency vectors ipk/lk, —1) which are orthogonal to the modes Uk = {qk,Pk) 
{uk is called "resonance" and is also denoted Pk'-<lk in what follows) accumulate at the invariant torus with frequency 
vector ijj — (cj,— 1). This family of periodic orbits satisfies the following relations : \uj ■ Uk+i] < \uj ■ Vk\ and 
limfe_,oo |a; • =0, and Vk is given by the following equation : 

Uk=Na,---Na,_,U0, (3.3) 

where i/q ~ (1,0) and Na^ denotes the matrix 



ai 1 
1 



The set of the Fourier modes of the perturbation with frequency vectors i>k leads to the divergence of perturbation 
expansions since the small denominators which are equal to a; • i/^ tend to zero as k increases. The renormahzation 
transformation deals with the modes specifically by non-perturbative techniques. 
The transformation TZ is composed by four steps : 

(1) A shift of resonances constructed from the condition Ui ^ : we impose that cos[(ao, 1) • v] = cos[(l,0) • 
if'] where ao = [^~^] is the integer part of uj^^ . This change of coordinates is performed by a linear canonical 
transformation 

{A,ip)^{A\^') = {N-^'A,Na,v), 

which is generated by F{A\ cp) = Na„A' ■ (p. The Hamiltonian expressed in the new coordinates becomes 

H\A', if') = H{A, ip)=ivA + V{n ■ A, cp), 

= iv ■ Na„A' + V{^1 ■ Na„A', iV,- V), 
= Na,u, ■ A' + V{Na,n ■ A', iV->')- 

Thus the new frequency vector is equal to Na^u} — —louj' where ui' — [uj' , —1) and uj' is the image of uj by the Gauss 
map 

UJ ^ uj' = uj-^ - [uj-\ (3.4) 
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This map corresponds to a shift to the left of the entries in the continued fraction expansion of the frequency 



uj = [ao,ai,a2, .. lo' ^ [ai,a2,a3, . . .]. 

The main effect of Step 1 is to change the frequency vectors of the Fourier modes of the perturbation according to 
the map v ^ N~^v. 

(2) We rescale the energy by a factor w"^ (or equivalently time by a factor uj), i.e. we multiply the Hamiltonian 
by w"^, and we change the sign of both phase space coordinates (A, tp) ^ (—A, —^), in order to have a;' as the new 
frequency vector, i.e. such that the linear term in the actions Nagi^ ■ A = —uju)' ■ A is rescaled into a term of the form 
uj' ■ A. Moreover, fl — (1, a) is changed into Nagft — (oq + a, 1) = (oq + a)ft' with ft' = (1, a') = (1, {qq + 
since the normalization conditio n we use is that the first component of n is equal to one. The map a t—f (oq + 



is the inverse of the Gauss map (3.4), in the sense that if a = [feo, ^i, • • ■ ] then a' — [ag, Bq, bi, . . . ]. If we define 

[a\uj] = [... ,62,6i,6o|ao,ai,a2, . . .], 

The map [a\uj] > [a'l'^'] corresponds to a two-sided Bernoulli shift : 

[• ■ • ,b2,bi,bo\ao,ai,a2, ■ ■ ■] ^ [■ ■ ■ ,b2,bi,bo,ao\ai,a2, ■ ■ ■]■ 

Since w < 1, the effect of the rescaling of time is that the trajectories of the rescaled Hamiltonian correspond to the 
ones of the initial Hamiltonian on a longer time scale. 

(3) We perform a rescaling of the actions : H is changed into 

H\A,v) = XH (^j,cp 

with A = X{H) such that the mean value of the quadratic term of H' in the variable fl' ■ A, is equal to {CI' ■ A)'^ /2 : 

X^2u;-\ao + a)^{V^^'>), (3.5) 



where V^^'^ denotes the coefhcient of the quadratic term of V{fl ■ A, ip) = H — liJ ■ A in the variabl e fl ■ A, and {V^ 
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denotes its mean value on [0, 27r[^. Since in general (F^^-') is close to 1/2, the rescaling coefficient (3.5) is larger than 
one. Thus the rescaling in the actions corresponds to a focus in phase space around the invariant torus with frequency 
u)' located approximately at fl' • A w 0. 
After these three steps, H is changed into 

H'iA, p)=.'.A + 2.-(ao + anvi^^)V (~ j ^)^^(2)^ ' ^> "^.oV 

(4) The fourth step is a canonical transformation that eliminates the non-resonant part (denoted /^) of the 
perturbation of H' . 

The choice of the part of the perturbation which has to be considered resonant or non-resonant is somewhat arbitrary. 
The set of non-resonant modes contains the modes of the perturbation which are sufficiently far from the resonances 
{I'k} in order to avoid small denominator problems during the elimination process. A convenient choice concerning 
the non- resonant modes is the set /~ of integer vectors v = (i^i, 1^2) such that |j/2| > \vi\ '■ 

r ={V={V^,V2)^1?\\V2\>\V1W. 



We notice that Eq. (3.3) defining i/^ — (qk,Pk) shows that pk < Qk for fc > 0. Consequently, the resonances do not 
belong to /~. At each iteration of the transformation, the frequency vector of the considered torus changes (since we 
perform unimodular transformations). We have chosen a unique region /" such that it does not contain any of the 
resonance lines u; • = for all cj e]0, 1[. 

From the form of the eigenvectors of Na^, we can see that each vector 1/ £ \ {0} is mapped into after a sufficient 
number of iterations of the matrices iVa. (the eigenvector of associated with the eigenvalue of modulus larger 
than one points into /^). In other terms, each resonant mode becomes non- resonant at a sufficiently smaller scale in 
phase space. We notice that is not an element of /~, i.e. it is resonant. 

Since the initial Hamiltonian (^]^) is quadratic in the actions, the renormalization for quadratic Hamiltonians defined 
in Refs. |^,^ is well-suited for this problem. We define more precisely Step 4 for the following Hamiltonians : 

H{A, p)=uj-A + m((p)(0 • A)^ + g{ip)fl ■ A + f{p), 
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where m, g and / are scalar functions of the angles and (m) ^ 0. We eliminate completely the non-resonant 
modes of g and / by a canonical transformation connected to the identity, which is defined by iterating KAM-type 
transformations. The KAM iterations we perform (by Lie transformations) are generated by functions that are linear 
in the actions. The Hamiltonian expressed in the new coordinates is again quadratic in the actions. Thus this type of 
transformation allows us to remain quadratic at each step of the transformation . One iteration Us of the KAM 
transformation reduces the non-resonant modes of / and g, from order e to e^. The transformation that eliminates 
completely the non-resonant part is defined in the following way : 

H' = H o Uh , where Uh = ° ° ' ' ' ^S„ o • • ■ , 

where the purpose of Us„ is to reduce the non- resonant part of / and g from order e^" ^ to e^" , such that /' = g' = 
0, where 1~ f denotes the non- resonant part of the constant term in the actions of H' , i.e. I^/' = X^i^e/- H^^'^ ■ 
For the explicit equations of this part of the transformation, we refer to Refs. fl^. 

In summary, the renormalization transformation acts in the following way : First, some of the resonant modes are 
turne d n on-resonant by a rescaling of phase space that changes the frequency of the torus according to the Gauss 



map (3.4). Then an iteration of a KAM-type transformation eliminates the non-resonant modes (by slightly changing 
the resonant ones). 

The renormalization transformation is a map of Fourier coefficients : If 



then H' = TZ{Ti.) becomes t/j = ^ {{Hk,i>}) where fc = 0, 1, 2 in our case. The properties of TZ are obtained by 
iterating this map. In order to analyze it numerically, we truncate this map, i.e. we truncate the Fourier series by 
considering frequency vectors v = {vi, V2) such that maxi=i.2 < L. We observe the convergence of the properties 
of the renormalization transformation as L grows. 



IV. CRITICAL FUNCTION OF THE MODEL 



A. Critical function e'c{u}) 

The aim is to obtain the value s'^{ijj) for which the invariant torus with frequency lo exists for e' < s'^^lu) and is 
broken for larger values. 

If Lu satisfies a Diophantine condition, the KAM theorem ensures the persistence of an invariant torus with 
frequency lo for the system we consider, i.e. s'^{uj) > 0. If w is rational, a resonance breaks up the torus, i.e. e'^{uj) = 0. 
Moreover, this function is symmetric : e^(— w) = s'^{ijj). This symmetry comes from the fact that the c anon ical 
transformation {Ai, ipi, A2, (^2) ^ (—^1, — ¥'1, ^2, ^2) only changes the frequency lo into —lo in Hamiltonian ( |2.3| ). 
The critical function e'^{iio) is determined by looking at the iterates of the renormalization transformation TZ described 
above, i.e. it is defined by the following equations : 

7^"i^e' ^ HQ{A)=^i:j-A+\[VL-Af fore'<4(tj), (4.1) 

n — >oc 2 

7^"i^e' ^ 00 for e' > 4(w), (4.2) 

n — >oo 



where H^i is Hamiltonian (3.1). Numerical work shows that Eqs.(4.1 )-(4.2) define uniquely £c(^); i-^- the renormal- 
ization has two main domains : one where the iterations of the transformation converge to an integrable Hamiltonian 
Hq — this domain is conjectured to be the set of Hamiltonians that have a smooth invariant torus with frequency 
u) — and a domain where the iterations diverge, and this domain is conjectured to be the set of Hamiltonians that 
do not have an invariant torus with the frequency uj. For e' = £c(^)i it was shown numerically that the iterations 
converge to a strange chaotic attractor containing all the relevant information on critical tori (i.e. at the threshold of 
the break-up) This means that we conjecture that the critical thresholds obtained by the complete renormalization 
coincide (up to numerical precision) with the thresholds of the break-up of invariant tori. In order to give support to 
this conjecture, we compare the values we obtain with another independent method, Greene's residue criterion. This 
method consists in analyzing the stability of nearby periodic orbits (for instance, the ones given by the truncations 
of the continued fraction expansion of the frequency of the torus). For rigorous results on Greene's criterion, we refer 
to Refs. Hfi. 
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Figure 1^ shows the value of e'^{uj) for to €]0, 1 [ determined by the renormahzation method with L = 20. A first remark 
is that the last invariant torus is not the golden mean one as it is the case for the standard map and for Escande's 
paradigm Hamiltonian : The torus with frequency 7 — (\/5 — l)/2 = [1, 1, 1, . . .] is broken for £^(7) ~ 0.02995 
whereas for instance the torus with frequency uji = (7 + 3)/5 = [1, 2, 1, 1, 1, . . .] persists until e'dt^i) ~ 0.03163. These 
two values coincide up to numerical precision with the ones obtained by Greene's residue criterion : Ec (7) ~ 0.0299 
and ef\u;i) « 0.0316. 

The last KAM torus is broken for e* w 0.03334. Its frequency is equal to UJ2 ~ 0.6976. The critical threshold of the 
break-up of this torus with frequency UJ2 obtained by Greene's residue criterion coincides with e* up to numerical 
precision : Sc (^2) « 0.0333. For e > e* , there are no invariant tori left, and large scale stochasticity occurs : 
trajectories can go from the resonance 1:1 (located at Ai = 1) to the resonance 0:1 (located at Ai = 0). 
Between two neighboring main resonances l:n and l:n-|-2, quasiperiodic motion with frequency uj €]l/(n + 2), l/n[ 
can occur for small e'. If e' is greater than some value e„, there is no longer any quasiperiodic motion in between 
these two resonances (and some chaotic trajectories can go from one resonance to the other). For n — 1, this value 
is equal to e* for which all rotational invariant tori are broken since the last invariant torus to break-up is located 



between the resonances 1:1 and 1:3. We apply Chirikov's criterion in order to have an estimate of £„ [see Eq.(2.5)] 

1 



£(^) 



4(n+ 1)2 



For the largest one n = 1, this value is equal to 0.0625 which is approximately twice the value obtained by renormal- 
ization. Escande's approximate renormalization gives 0.0352 as the critical amplitude of the field for the last invariant 
torus, which is close to the value determined by the complete renormalization method. This feature is expected to be 
true for quadratic Hamiltonians in the actions as it was pointed out in Refs. ||20| , pT|] . In Ref. [^, it has been noticed 
that the approximate renormalization usually slightly overestimates the re al c ritical value. 



Lin and Reichl |ll[| have developed a method adapted to the specific model (IT) to compute the critical amplitudes e„. 
This method is based on the fact that as soon as the last invariant torus is broken, some trajectories starting near one 
of the main resonance (say l:n) can approach the other main resonance l:n-|-2. Since the diffusion of these trajectories 
can be very slow (due in particular to the resonances of low order between two main neighboring resonances, and 
in particular in the region of phase space where the last KAM torus breaks up when the parameter e' is close to its 
critical value), the values they obtained overestimate the ones obtained by renormalization. For instance, between 1:1 
and 1:3 the critical value they obtained is approximately 0.037±0.001. 

(c) 

For n ~ 3, the estimate obtained by Chirikov's criterion is £3 « 0.156, Escande-Doveil's renormalization gives 
£3 « 0.0080, whereas Lin and Reichl obtained £3 « 0.0081 ± 0.0003. The approximate value of £3 obtained by the 
complete renormalization is £3 « 0.0068. Again, the same comments apply to this case : Lin and Reichl's method 
overestimates the value obtained by renormalization because of the very slow diffusion of the trajectories between the 
resonances 1:3 and 1:5 for critical or near-critical values of the parameter e' . 



B. Critical function ec{u};m, a,Q) 

We have studied the critical function e'^{uj) of Hamiltonian ( |2.3D . This rescaled Hamiltonian is equivalent to the 
initial Hamiltonian ( |2.2| ) with maO,^ = 1 and in particular with a = 1, m = 1 and f2 = 1. For other values of a, m, 
and Q, the critical function £c(w; m, a, fl) is equal to 

£c(w; m, a, f2) = maVt^e'^ , 

according to Eq. ( p.4[ ). We notice that the argument of e'^ is the rescaled frequency uj/Vt since time has been rescaled 
by a factor 17. Thus the critical function varies like the square of the frequency of the field. In particular, the largest 
value of the parameter for which an invariant torus persists between the resonances 1:?! and l:n+2 varies like : 

en[n) = n\n{n = 1), 

which is consistent with the numerical results found in Ref. ijll]] . This feature can be generalized to the following 
one-parameter family of Hamiltonians : 

2 
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where 5, are 27r-periodic functions. This can be seen by rescahng time by a factor CI 



and the actions by a factor 1 /SI 



H"{p',x,t') = ^H\np',x,t'). 



The fact that the critical couphng Sc{i-o: m, a, CL) is proportional to m is general for a particle in an infinite square- well 
potential driven by a perturbation depending on x and periodically on t. This is obtained by rescaling the momentum 
by a factor 1 /m of the Hamiltonian : 

2 

H{p,x,t) = ^+VsQ{x)+eV{x,t). 

From the fact that the interaction with the field is proportional to x (the interaction is of the form ex f (fit) where / 
is 27r-periodic), the critical function Sc{co;m,a,Q) is expected to be proportional to a. 



V. CONCLUSION 



Wc have defined and studied numerically a rcnormalization transformation for a system of a particle in an infinite 
square well potential driven by an external monochromatic field in order to determine the critical thresholds of the 
break-up of invariant tori. Rcnormalization allows to obtain very precise information on the stability of Hamiltonian 
systems with two degrees of freedom, and in particular it allows to determine what kind of stable motions remains as 
a function of the amplitude of the perturbation. We have checked that for some specific frequencies, the thresholds 
obtained by rcnormalization coincide with the ones obtained by Greene's residue criterion. 

We have chosen to study a particle in an infinite square well potential driven by an external field because of its 
simplicity, but the renormalization methods are very general and can be applied directly to other types of Hamiltonian 
systems with two effective degrees of freedom. 
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FIG. 1. Critical function of Hamiltonian system (B.l) 
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